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Experimental System — CVB Experiment 
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• Partially filled cell forms a Constrained Vapor Bubble (CVB) design. 

• Wickless heat pipe with transparent walls. 

• Goal is to measure the vapor and liquid distributions as the pipe is driven. 
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Evaporation from a Thin Film 
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Transition region controls evaporation. Lowest overall resistance. 

- Accurate modeling will enable us to engineer optimal surfaces 
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Reflectivity / Interferometry T echnique 



Varying thickness of the meniscus 
produces an interference pattern 


Analysis yields: 

Film thickness profile 
Contact angle 
Interface curvature 
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Interference pattern analyzed to obtain 
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Evaporating Meniscus — HFE-7000 
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Fluid Flow Model 




Lubrication approximation used to model fluid flow. 


Navier slip (solid-liquid interface) and Marangoni 
shear (liquid-vapor interface) boundary conditions 
applied 
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• Temperature dependence of fluid properties accounts for the capillary, Marangoni and van 
der Waals forces. 

• Mass balance, provides the evaporating mass flux at each pixel location. 
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Heat Transfer at the Contact Line 



• Heat transfer at the contact line was modeled using a Kelvin-Clapeyron approach. 
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• The temperature difference across the liquid was tied to the mass flux to describe 

the driving force as a superheat. 
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• Final equation can be written as a 4th order differential equation for the film 
thickness. 
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Comsol Set-Up 



• Split the 4 th order equation into two 
2 nd order equations for film 
thickness and interface curvature. 

• 1-D weak form on the boundary. 

• 2-D conduction heat transfer in the 
solid. 

• Boundary conditions set film 
thickness and curvature at the thick 
end, and the curvature and slope at 
the start (perfectly wetting fluid). 
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Simulation Results — Steady-State 



COMSOL simulation was able to 
reproduce the experimental data 
from an octane meniscus. 
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hydrodynamic slip was necessary « 

to match both the position of the | 

curvature peak and the spread of j= 

the peak. i§ 

Peak height is controlled by the 
thickness of the adsorbed film 
ahead of the contact line. 
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Simulation Results — Oscillation 



• COMSOL simulation was applied to an 
pentane meniscus. 

• Oscillation was established by varying 
the temperature on the underside of the 
substrate. 

• Plots show snapshots of the temperature 
profile in the substrate for two different 
substrate materials at the same time 
during the oscillation. 

• Higher thermal diffusivity leads to a 
different dynamic mismatch between 
heat delivered to the film and heat 
removed via evaporation. 
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Dimensionless Heat Flow Rate 


Simulation Results — Oscillation 



Achieve large changes in film thickness 
for small differences in applied 
temperature (here -0.01 °C) 




Film thickness at evaporation front 



• Dimensionless heat flow rate changes 
dramatically with substrate material. 
The oscillation phase also changes 
and for certain forcings, the film 
oscillates at higher frequency. 
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Simulation Results — Oscillation 
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• Film acting as a low-pass filter. 

• There is a small resonance that we need to explore further. 
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Conclusions 



■ We can reproduce some of what we observe in an oscillating 
evaporating meniscus. 


■ Oscillation amplitude and frequency are related to a mismatch 
between the heat dissipation rate in the film and the rate at 
which heat can be delivered from the substrate. 


■ Much more modeling work is needed to match up with actual 
experimental data and to mimic all the behaviors we observe in 
the thin film. 
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